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Abstract 

We present numerical calculation of the evolution of a background space-time metric and sub- 
horizon matter density perturbations in viable f{R) gravity models of present dark energy and 
cosmic acceleration. We found that viable models generically exhibit recent crossing of the phantom 
boundary u?de = — 1- Moreover, as a consequence of the anomalous growth of density perturbations 
during the end of the matter-dominated stage, their growth index evolves non-monotonically with 
time and may even become negative. 
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I. INTRODUCTION 



It is one of the most important issues for cosmologists and particle physicists to under- 
stand the physical origin of the dark energy (DE) which is responsible for an accelerated 
expansion of the current Universe. Although the standard spatially flat A-Cold-Dark-Matter 
(ACDM) model is consistent with all kinds of current observational data 1], some tentative 

nn \ 

deviations from it have been reported recently |2|, |3J. Furthermore, in the ACDM model, 
the cosmological term is regarded as a new fundamental constant whose observed value is 
much smaller than any other energy scale known in physics. Hence it is natural to seek for 
non-stationary models of the current DE. Among them, f(R) gravity which modifies and 
generalizes the Einstein gravity by incorporating a new phenomenological function of the 



ar R, f(R), can provide a self-consistent and non-trivial alternative to the ACDM 
5L 



Ricci sea 
model 

In the previous paper [6J, we calculated evolution of matter density fluctuations in viable 
f(R) models [4, 5] for redshifts z 3> 1 during the matter-dominated stage and found an 
analytic expression for them. In this paper we extend the previous analysis and perform 
numerical calculations of the evolution of both background space-time and density fluctu- 
ations for the particular f(R) model of Ref. ^] without such a restriction. As a result, we 
have found crossing of the phantom boundary u^e = —1 at an intermediate redshift z < 1 
for the background space-time metric and an anomalous behavior of the growth index of 
fluctuations. 



II. BACKGROUND 



We adopt the following action of f(R) models with model parameters n, A and R 6 



S = TF~n I d A x^f(R) + S m , f(R) = R + \R S 
IbnG J 
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(i) 



where S m is the action of the matter content which is assumed to be minimally coupled to 
gravity. To make the late-time asymptotic de Sitter regime where R = constant stable, A 
has to satisfy f'(R) > Rf"(R). As a result, A has a lower limit A m ; n for each n. Numerically 
we find (n, A min ) =(2, 0.9440), (3, 0.7259), and (4, 0.6081). From the action flQ, we obtain 
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field equations as 

K - \KR = -SvrG (Vj m) + T; (DE) ) , (2) 
SttGT^ = (F- 1)J2JJ - i(/ - fl)^ + (VV„ - #D)F, F(J2) = f'(R). (3) 

Working in the spatially flat Friedmann-Robertson- Walker (FRW) space-time with a scale 
factor a(t), 

3H 2 = 8nGp - 3(F - l)H 2 + -(FR - f) - 3#F, (4) 
2H = -8nGp - 2(F — l)H — F + HF, (5) 

where H is the Hubble parameter and p is the energy density of the material content which 
we assume to consist of non-relativistic matter. From fll]), we can determine the DE equation 
of state parameter % E , 

Pde 1 ^ 2H(F-1)-HF + F 

W DE = = — H : : • (6) 

Pde -3HRF' + 3(H 2 + H)(F-l)-(f-R)/2 KJ 

We solve the evolution equation (jSJ) numerically using (jl]) to check numerical accuracy. 
The moment when the matter density parameter Q(t) = lQirGp/ (IQnGp + XR S ) equals 
to 0.998 is chosen as the initial time £j. We determine the current epoch t = to by the 
requirement that the value of Q takes the observed central value Qq = 0.27. R s is fixed in 
such a way as to reproduce the current Hubble parameter Hq = 72km/s/Mpc. We find the 
ratio R s /Hq is well fit by a simple power-law R s /Hq = c n X~ Pn with (n,c n ,p n ) =(2, 4.16, 
0.953), (3, 4.12, 0.837), and (4, 4.74, 0.702), respectively, whereas in the ACDM limit it 
would behave as R s /H 2 = 6(1 - fi )/A ~ 4.38A- 1 . 

Figures [1] depict evolution of u>de as a function of redshift z where phantom crossing is 
manifest. As expected, it approaches wde = — 1 = constant as we increase A for fixed n. 
For the minimal allowed values of A, deviations from wde = — 1 are observed at 5% level on 
both directions in z < 2 independently of n. From ([6]), we can read off that this phantom 
crossing behavior is not peculiar to the specific choice of the function ([T|) but a generic one 
for models which satisfy the stability condition F' > 0. 
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FIG. 1: Evolution of the equation-of-state parameter of effective dark energy. 
III. PERTURBATIONS 

We now turn to the evolution of density fluctuations. In f(R) gravity, the evolution 
equation of density fluctuations, 5, deeply in the sub-horizon regime is given byp| 

'5 + 2H5-4<kG cSP 5 = 0, G cS = - (7) 

1 + 6 a?T 

In the previous paper jg] we obtained an analytic solution for the high-curvature regime 
when the scale factor evolves as a(t) oc t 2 ^ 3 and F takes an asymptotic form F ~ 1 — 
2nA (R/ R s )~ 2n ~ . In the present paper, we numerically integrate (j7]) up to z — without 
using the approximation \F — 1| <C 1. 

The wavenumber of our particular interest is the scale corresponding to ag normalization, 
for which we find k e g(r = 8/i -1 Mpc) = 0.174/iMpc -1 . Since the standard ACDM model 
normalized by CMB data explains galaxy clustering at small scales well, 5fRG should not be 
too much larger than #acdm on these scales. We may typically require (5fRc/^ACDM) 2 ^ l l- 
Although we neglect non-linear effects here, the difference between linear calculation and 
non-linear N-body simulation remained smaller than 5% at wavenumber 0.174/iMpc -1 8]. 

The left panel of Fig. [2] represents ((SfRG/^ACDM) 2 as a function of A for n = 2 together 
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with two fitting functions. The solid line is from the analytic formula obtained in Ref. 
and the broken line is numerical fitting using an exponential function 1 + b n e~ qnX . From 
these analysis, we can constrain the parameter space as the right panel of Fig. |2j The region 
which satisfy (<5fRc/ <5acdm) 2 < 1.1 lies above the solid line. The region below the dotted line 
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FIG. 2: Constraints from the power spectrum. 



is forbidden from instability of de Sitter regime. 

Next we turn to another important quantity used to distinguish different theories of 
gravity, namely, the gravitational growth index, 7(2), of density fluctuations. It is defined 
through 
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n m ( z y^, or 7 (z) 



l°g life 



(8) 



(I In a \ogfl m 
It takes a practically constant value 7 = 0.55 in the standard ACDM model, but it evolves 
with time in modified gravity theories in general. We also note that 7(2;) has a nontrivial 
/c-dependence in f(R) gravity since density fluctuations with different wavenumbers evolve 
differently. Therefore, this quantity is a useful measure to distinguish modified gravity from 
ACDM model in the Einstein gravity. 

Figures [3] show the evolution of 7(2;) together with that of G e ^/G for different values of k. 
In the early high-redshift regime, 7(2) takes a constant value identical to the ACDM model. 
It gradually decreases in time, reaches a minimum which may be even negative, and then 
increase again towards the present epoch. We note that recently Narikawa and Yamamotoj^] 
numerically calculated time evolution of 7(2) in a simplified model which we had used in 
the previous paper and also obtained some analytic expansion, which behaves qualitatively 
the same as the numerical result but with much more exaggerated amplitudes. Our results, 
which satisfy all viability conditions, exhibit milder deviation from ACDM model than they 
found. At present, the constraints for the growth index is not so strict to distinguish the 
deviation from the ACDM model lOj, but observations may reveal its time and wave number 
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FIG. 3: Evolution of and G cS /G. 

dependence in future. 

IV. CONCLUSION 

In the present paper we have numerically calculated the evolution of both homogeneous 
background and density fluctuations in a viable f(R) model of accelerated expansion based 
on the specific functional form proposed in Ref. [5j]. We have found that viable f(R) gravity 
models of accelerated expansion generically exhibit phantom behavior during the matter- 
dominated stage with crossing of the phantom boundary u>de = — 1 at redshifts z < 1. The 
predicted time evolution of u>de has qualitatively the same behaviour as that was recently 
obtained from observational data in [2j . 

As for density fluctuations, we have numerically confirmed our previous analytic results 
on a shift in the power spectrum index for large wavenumbers which exceed the scalaron 
mass during the matter dominated epoch[6|, while for smaller wavenumbers fluctuations 
have the same amplitude as in the ACDM model. 

We have also investigated the growth index j(k, z) of density fluctuations and have given 
an explanation of its anomalous evolution in terms of time dependence of G e Q. Since 
7 has characteristic time and wavenumber dependence, future detailed observations may 
yield useful information on the validity of f(R) gravity through this quantity, although 
current constraints have been obtained assuming that it is constant both in time and in 
wavenumber (2!, [lo| . 



6 



Acknowledgments 



HM and JY are grateful to T. Narikawa and K. Yamamoto for useful communications. AS 
acknowledges RESCEU hospitality as a visiting professor. He was also partially supported by 
the grant RFBR 08-02-00923 and by the Scientific Programme "Astronomy" of the Russian 
Academy of Sciences. This work was supported in part by JSPS Grant-in-Aid for Scientific 
Research No. 19340054(JY), JSPS Core-to-Core program "International Research Network 
on Dark Energy", and Global COE Program "the Physical Sciences Frontier", MEXT, 
Japan. 



[1] E. Komatsu et al. larXiv:1001.4538l 

[2] A. Shafieloo, V. Salmi and A. A. Starobinsky, Phys. Rev. D 80, 101301 (R) (2009) 

|arXiv:0903.5141j . 
[3] R. Bean, arXiv:0909.38531 

[4] W. Hu and I. Sawicki, Phys. Rev. D 76, 064004 (2007) [arXiv:0705.fl58] . 
[5] A. A. Starobinsky, JETP Lett. 86, 157 (2007) [i arXiv:0706.2041| . 

[6] H. Motohashi, A. A. Starobinsky and J. Yokoyama, Int. J. Mod. Phys. D 18, 1731 (2009) 

|arXiv:0905.0730| . 
[7] S. Tsujikawa, Phys. Rev. D 76, 023514 (2007) [arXiv: 0705 .1032] . 
[8] H. Oyaizu, M. Lima and W. Hu, Phys. Rev. D 78, 123524 (2008) [ arXiv:0807.2462] . 
[9] T. Narikawa and K. Yamamoto, larXiv:0912. 14451 
[10] D. Rapetti, S. W. Allen, A. Mantz and H. Ebeling. la7Xiv:0911.1787l 



7 



